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heterostructures: A density functional theory study 
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A computational study of tiie epitaxial Co2MnSi(001)/MgO(001) interface relevant to tunneling 
magnetoresistive (TMR) devices is presented. Employing ab initio atomistic thermodynamics, we 
show that the Co- or MnSi-planes of bulk-terminated Co2MnSi form stable interfaces, while pure Si 
or pure Mn termination requires non-equilibrium conditions. Except for the pure Mn interface, the 
half-metallic property of bulk Co2MnSi is disrupted by interface bands. Even so, at homogeneous 
Mn or Co interfaces these bands contribute little to the minority-spin conductance through an MgO 
barrier, and hence such terminations could perform strongly in TMR devices. 



In recent years, tunneling magneto-resistance (TMR) 
devices, consisting of two ferromagnetic electrodes sepa- 
rated by a tunneling barrier in the form of a thin oxide 
layer, have been subject to intense research. This is due 
to the use of TMRs in non- volatile magnetic random ac- 
cess memories (MRAMs), which are expected to play an 
increasing role in future electronics. The TMR device 
allows for electric read-out of magnetically stored infor- 
mation, as the tunneling resistance is strongly dependent 
on the relative orientation of the magnetization in both 
electrodes (parallel or anti-parallel). The difference be- 
tween the electrical conductances for these two orienta- 
tions, divided by the smaller of the two, called the TMR 
ratio, could possibly be huge if the material used for one 
or both electrodes is a ferromagnetic half-metal. In such 
a material, the Fermi energy Ep lies in a gap for the 
electrons of one spin, while the density of states (DOS) 
is metallic for the electrons of the other spin. Examples 
of present interest are the Heusler alloys Co2MnSi and 
Co2FeSi (and variants thereof), materials that are pre- 
dicted by calculations P, 0,01 to be bulk half- metals. In a 
real TMR device, spin-flip processes related to electronic 
interface states could invalidate the advantages of half- 
metallic electrodes. [41] Even so, spin-selective tunneling 
may help to operate these devices with high TMR ratio. 

Experimentally, progress towards devices with large 
TMR ratio has been achieved by using MgO tunnel- 
ing barriers. The success of this barrier material was 
traced back to selection rules based on band symme- 
try that may greatly enhance the ratio of the spin- 
polarized currents. 0, Q By combining an MgO barrier 
with Co2MnSi electrodes, TMR ratios as high as 192% 
and even 753% [1] have been achieved. However, the 
underlying materials issues limiting the achievable TMR 
ratio have largely remained elusive, mostly due to the 
unknown structure of the ferromagnet-oxide interface. 

Here, we employ density functional theory (DFT) cal- 
culations to provide comprehensive information about 
the structural stability, the magnetic and the elec- 



tronic properties of the Co2MnSi(001)/MgO(001) inter- 
face. The highly accurate all-electron full-potential lin- 
earized augmented plane wave (FP-LAPW) method im- 
plemented in the Wien2k code is used to calculate 
total energies and electronic structure of the interfaces. 
The generalized gradient approximation (PBE96 [l^) to 
the exchange-correlation functional has been adopted be- 
cause the GGA gives better results than the LDA (re- 
garding lattice constants and bulk moduli) for Co2MnSi 
For bulk Co2MnSi , the GGA-PBE calculations 



give us a gap in the minority-spin DOS of 0.84eV, and 
a spin gap (between the top of the minority valence 
band and Ep) of 0.33eV [l2j . in excellent agreement 
with previous calculations [3[. The muffin-tin radii were 
set to i?MT(Mg) = Rmt{0) = 1.8 bohr, i?Mr(Co) = 
i?Mr(Si) = 2.0 bohr and i?Mr(Mn) = 2.1 bohr. A k- 
point mesh of 16 x 16 x 2 and in the interstitial region 
a plane wave expansion cutoff of 21 Ry was used. All 
internal atomic coordinates were relaxed until the forces 
on the nuclei were smaller than 2 mRy/bohr. 

The TMR elements are modeled in supercell geome- 
try, using 3 layers of MgO and 11 or 13 layers of the 
Heusler alloy, depending on the termination at the inter- 
face. The supercell is chosen symmetric with respect to 
a mirror plane at the central layer of the MgO barrier, 
with two identical interfaces. Epitaxial growth of the 
heterostructure is assumed. We use a (1 x 1) unit cell 
with the calculated lattice constant of Co2MnSi, 5.629A 
parallel to the interface, and allow for a tetragonal dis- 
tortion of the thin MgO film. (The bulk lattice constants 
of Co2MnSi and MgO differ by 7% in PBE-GGA). The 
length of the supercell in (001) direction is optimized for 
each interface in a series of calculations, thus allowing for 
the bond distance at the interface to adjust. 

First, we study the thermodynamic stability of the 
Co2MnSi/MgO interface. Four different terminations of 
Co2MnSi(001) (by CoCo, MnSi, MnMn or SiSi planes), 
each in combination with four different registries rela- 
tive to MgO(OOl) (O top, Mg top, bridge, or hollow site) 
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FIG. 1: (Color online) Left: Ball-and-stick models of the re- 
laxed atomic structures of the MnMn/0 interface (top), the 
CoCo/O interface (middle), and the MnSi/0 interface (bot- 
tom). The numbers give the spin magnetic moment of the 
atoms in units of the Bohr magneton. Right: Phase diagram 
for the epitaxial Co2MnSi(001)/MgO(001) interface. The col- 
ored regions correspond to different interface terminations be- 
ing stable under the conditions described by the chemical po- 
tentials Afii = fii — Qi, where gi {i = Co, Mn) is the total 
energy of one unit cell of the stable bulk phase of the ele- 
ments (hep Co and a-Mn). The polygon bounded by the red 
dashed lines indicates the region accessible in thermodynamic 
equilibrium with the bulk phases of Mn, Co, MuaSi or CoaSi. 



have been taken into consideration. By comparing the 
DFT total energies, we find that for the CoCo, MnSi 
and MnMn terminated Co2MnSi the O top site is lower 
in energy than any other registry, and thus stable. For 
these three cases, the Mg top site is metastable, while 
the bridge site is unstable, and the atoms relax to the O 
top site instead. For the SiSi termination, the Mg top 
site is found to be the most stable position. 

In order to compare the stability of different termina- 
tions, we use the method of ab initio atomistic thermo- 
H [Hi. The interface energy, normalized to 



dynamics 
the area A 



AGiN,,T,p)-Y,^N,fi,{T,p) 



(1) 



Here, AG is the difference in Gibbs free energy between 
the heterostructure and appropriate amounts of the two 
bulk materials forming the interface. AN^ {i = Co, Mn, 
Si) counts the number of atoms in which the supercell 
deviates from bulk composition, and Hi is the chemical 
potential of element i. Since the interface energy is a 
free-energy difference, it is often a good approximation to 
calculate 7 from the differences of total energies obtained 
from DFT calculations. Differences in the other contri- 
butions to G (due to, e.g., the free energy of lattice vibra- 
tions) largely cancel, and are therefore neglected For 
the thermodynamic analysis, it is assumed that the in- 
terfaces are in thermodynamic equilibrium with the bulk 
Heusler alloy. This allows to eliminate one chemical po- 
tential via the relation 2/ico + A^Mn + A^Si = ffCosMnSi- 
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TABLE I: Spin polarization Pint of the density of states at 
the Fermi level, calculated by projecting onto the orbitals of 
the interface and sublayer atoms, for four terminations of the 
Co2MnSi/MgO(100) interface. Evbo is the band offset be- 
tween the valence band of MgO and the minority-spin valence 
band of Co2MnSi. 



Here, gcoaMnSi is the Gibbs free energy of Co2MnSi per 
formula unit. The interface energy thus depends on two 
independent chemical potentials, which we chose to be 
/iCo and A*Mn, as in previous studies. The interface phase 
diagram (Fig. [T|) displays the termination with the min- 
imal interface energy for every pair (/iMn,MCo)- 

We conclude from the phase diagram shown in Fig. [T] 
that both the termination of the Heusler electrode by a 
Co plane, sitting on top of oxygen, as well as termina- 
tion by a MnSi plane, with both species sitting on top of 
oxygen, are thcrmodynamically stable, depending on the 
preparation conditions of the Co2MnSi(001)/MgO(001) 
interface. The other two terminations considered are not 
accessible within the limits set by thermodynamic equi- 
librium with other phases. However, they could possibly 
be fabricated using out-of-equilibrium conditions, and 
therefore will be considered further. Our DFT calcula- 
tions show sizeable structural relaxations at the MnSi- 
intcrface (Fig. [1] lower left). An Mn-0 bond with a 
length of 2.27A is formed, while the Si atoms relax to- 
wards the Heusler alloy and shorten their bonds to the 
Co atoms in the next-nearest crystal plane (sublayer). 
In contrast, the Co- and Mn-terminated interfaces re- 
main planar, with Co-0 and Mn-0 bond length of 2.09A 
and 2.40A, respectively (Fig. [U middle and upper left). 
These results agree with previous DFT studies using a 
pseudopotential approach. [3| 

Magnetic moments of interface atoms are generally 
very sensitive to atomic structure. Our calculations show 
that only the spin magnetic moments at the interface and 
in the first sublayer of the Heusler alloy (cf. Fig. [T]) devi- 
ate significantly from the bulk moments, mco ~ 1.07/iB, 
TOMn — 2.91/iB, and mgi = — 0.04/1^. Comparison of 
the calculated moments in Fig. [l]to experimental values 
can be helpful in drawing conclusions about the interface 
chemistry. Measurements using X-ray magnetic circular 
dichroism on two Co2MnSi films of different thickness 
pointed to an increase of the magnetic moments at the 
interface to MgO. 17| Agreement with this experiment is 
closest if we assume a MnMn/0 interface (despite the 
Co- rich film composition reported in Ref. \17\ ). as only 
this structure shows enhanced Co moments. 

The interplay between the atomic and the electronic 
structure of the interfaces is best analyzed by inspecting 
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FIG. 2; (color online) The orbital-resolved DOS of the atoms 
at the interface (left column) and in the first sublayer (right 
column) of the Heusler alloy for the three non half-metallic 
heterostructures CoCo/O, MnSi/O and SiSi/Mg (rows). The 
second label in each sub-plot denotes the species that is pro- 
jected out. Red (full) lines show the projected DOS for the 
'out-of-plane' orbitals {d^^2_^2 + d^z + rfyz for Co and Mn, 
Pz for Si). Blue (dashed) lines show the projected DOS for 
the 'in-plane' orbitals {d^2_y2 + dxy and + Py)- The black 
(dotted) line is a projection onto the Si-s orbital. 



the orbital-resolved densities of states (DOS). In Fig. [21 
the orbital-resolved DOS in the minority spin channel is 
shown for the three non-half-metallic interfaces. It is ob- 
served that the CoCo/0 interface has only a small DOS 
at Ep which involves both in-plane and out-of-plane Co 
orbitals, as well as out-of-plane orbitals of the sublayer 
Mn atoms. The MnSi/O interface shows a sharp peak 
of the DOS at Ep due to an interface state. The or- 
bital analysis reveals that it arises from the out-of-plane 
orbitals of the sublayer Co atoms that hybridize with 
out-of-plane Mn orbitals in the interface layer, similar 
to the surface states studied previously [ll'|. In case of 
the SiSi/Mg interface, the half-metallic gap is completely 
filled with interface states originating from Si orbitals. 
These interface states may drastically reduce the spin 
polarization of the DOS at the Fermi level, Pint in Ta- 
ble m We find that the spin polarization breaks down 
(and even changes sign) for the MnSi/O and SiSi/Mg 
interfaces. Only the MnMn/0 interface is half-metallic, 
while for the CoCo/0 interface, a rather high degree of 
spin polarization, Pint = 67%, is retained despite the in- 
terface state. Similar to the Co2MnSi(001) SMr/ace[ll| , 



the strong hybridization of d^z and dyz orbitals of both 
first-layer Mn and second-layer Co atoms prevents the 
splitting-off of an interface state in this case. 



The results for the valence band offsets Eybo reported 
in Table [J are calculated from the self-consistent Kohn- 
Sham potential (see Ref. [3 for the approximations in- 
volved). The electronic core levels are used to define 
a common reference energy [l^. All interfaces show a 
sizable valence band offset, with the exception of the 
SiSi/Mg interface, where the two valence bands are al- 
most perfectly aligned. Hence, the SiSi/Mg interface is 
clearly unsuitable for TMR devices. Using as input i?vBO 
and the known band gap of MgO, we conclude that the 
band offsets for the conduction bands are > IV for all 
interfaces considered. This is different from the previous 
investigations of the interfaces of Co2MnGe with GaAs 
or Ge, where DFT calculations found Fermi-level pin- 
ning by gap states, or alignment of the valence bands, 
respectively [i^. 

Finally, we discuss the role of the interface states for 
the electrical properties. In order to achieve a large TMR 
ratio, it is crucial that the current through the MgO bar- 
rier is minimal in the state of anti-parallel magnetization 
of the two electrodes. For perfect half-metallic electrodes, 
this current results from a sequence of spin-flipping and 
tunneling to or from an interface state. We make the 
(conservative) assumption that the latter process is rate- 
limiting (cf. the discussion in Ref. [3)- To estimate the 
TMR ratio, one needs to take into account that the tun- 
neling probability of carriers in the various bands may 
be vastly different depending on band symmetry. For 
MgO, it has been demonstrated that band states of Ai 
symmetry show a weak decay inside the MgO barrier, 
while transmission of states with different symmetry is 
exponentially suppressed. 1, d For CoaMnSi, DFT cal- 
culations show that a conduction band in the majority 
spin channel with Ai symmetry exists, jl^ Hence, trans- 
mission of spin-majority electrons will be easy for paral- 
lel magnetization of the TMR electrodes, which is one of 
the pre-requisites for achieving a high TMR ratio. The 
crucial point, however, is the efhcient suppression of cur- 
rents for anti-parallel magnetization, which requires the 
tunneling rate into the minority interface states to be 
small. From our calculations, we can estimate this rate 
from the energetic splitting between bands having even 
and odd symmetry with respect to the mirror plane in- 
side the MgO barrier. In Fig. [3l the minority spin band- 
structures are shown. It is seen that the interface state 
at MnSi/0 shows considerable splitting at F, while the 
interface states at CoCo/0 and MnMn/0 show only a 
very small splitting of 41 meV and 11 meV, respectively 
(arrows in Fig. [3|). As seen from Fig. [21 the interface 
band at MnSi/0 consists of 3d orbitals perpendicular to 
the interface, which are invariant under symmetry op- 
erations in the interface plane, and thus match the Ai 
band of MgO; hence the large transmission. For the 
CoCo/0 and MnMn/0 interfaces, the interface bands 
formed by the in-plane 3(i orbitals transform under var- 
ious S-representations of the in-plane symmetry group; 
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FIG. 3; Minority spin band structures of the TMR elements. Only for the MnMn/O interface half-metallicity is retained. In 
the other cases interface states appear at the Fermi energy. The splitting of interface states with even and odd symmetry due 
to tunneling is indicated by the arrows. 



therefore the coupling to the Ai-type MgO band and the 
interface band splitting are small. As expected for a split- 
ting due to tunneling, it is largest at the F point, since 
the tunneling probability is largest for vanishing paral- 
lel crystal momentum. Knowledge of the band struc- 
ture puts us in position to speculate about the expected 
current-voltage characteristics of a TMR devices. For 
the MnMn/O interface, a bias voltage of ^ 0.2eV must 
be applied before electrons can be injected into the un- 
occupied interface band. For the CoCo/0 interface, the 
splitting of the interface state, and hence the tunneling 
rate of spin- minority electrons, is tiny at Ep. Only af- 
ter applying a bias voltage of ~ 0.35 eV, raising Ep in 
the spin-majority electrode, efficient tunneling into the 
spin-minority interface state near F will be energetically 
allowed. In TMR devices, the step-wise increase of the 
anti-parallel tunneling conductance as a function of bias 
voltage, as observed e.g. in Ref.d and 21, at low temper- 
atures, may be an indication of the role of interface states 
in transport. Then, the achievable TMR ratio in present 
TMR devices would still be limited by the (improvable) 
interface quality, rather than by factors intrinsic to the 
Heusler alloy, as suggested by Ref. [2l|. 

In summary, our DFT calculations indicate that TMR 
devices fabricated from Co2MnSi/MgO(100) heterostruc- 
tures may form thermodynamically stable CoCo/O or 
MnSi/0 interfaces, while the (most desirable) half- 
metallic MnMn/O interface may be prepared under non- 
equilibrium conditions. We conclude that interfaces 
made from a single atomic species (Co or Mn) enable 
a high TMR ratio, since the tunneling of the associated 
electronic interface states in the minority-spin channel 
through the MgO barrier is weak, which is a pre-requisite 
for a tolerably low parasitic current in the anti-parallel 
magnetization of the TMR device. 
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